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Abstract 

We study the steady state properties of a genotype selection model in presence 
of correlated Gaussian white noise. The effect of the noise on the genotype selection 
model is discussed. It is found that correlated noise can break the balance of gene 
selection and induce the phase transition which can makes us select one type gene 
haploid from a gene group. 
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1 Introduction 

Recently, nonlinear stochastic systems with noise terms have attracted extensive inves- 
tigations. The concept of noise-induced transition has got wide applications in the field 
of physics, chemistry and biology [1][2][3]. In most of these theories the noise affects the 
dynamics through system variable, i.e., the noise is multiplication in nature[4]. The focal 
theme of these investigations is steady state properties of systems in which fluctuations, 
generally applied from outside, are considered independent of the system's characteris- 
tic dissipation. These studies consider only one type of fluctuations that affect system's 
parameters, i.e., they drive multiplicatively. However, fluctuations due to some extrin- 
sic factors affect the system directly , i.e., they drive system dynamics additively. Since 
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two types of fluctuations have the common origin, they correlated with each other in the 
relevant timescale of the problem [5]. On the level of a Langevin-type description of a 
dynamical system, the presence of correlation between noise can change the dynamics of 
the system [6] [7]. Correlated noise processes have found applications in a broad range of 
studies such as steady state properties of a single mode laser [8], bistable kinetics [9], di- 
rected motion in spatially symmetric periodic potentials [10], stochastic resonance in linear 
systems[ll], and steady state entropy production[12]. In this paper we study the genotype 
selection model in presences of the correlated noise and discuss how noise correlation can 
break the balance of the gene selection and induce the phase transition. 



2 The simple genotype selection model 

We select a haploid group as our object and suppose that each haploid may have gene A 
or gene 5 [13]. The number of gene A haploid, gene B haploid and the total are N A ,N B 
and iV,respectively. Because the total number is constant (N = N A + A r ^=constant), we 
can make a transform as follow: 

x = —,l~x = —,0<x<l. (1) 

x is the ratio of gene A number to the total number. If gene A haploid and gene B 
haploid may have a mutation (A — ► B or B — ► A) during the process of the heredity 
and the ratio of A — ► B and B — ► A are m^At and mgAt, respectively, we can get the 
following difference equation: 

x{t - At) - x(t) = -m A Atx{t) + m B At[l - x(t)] (2) 
On the other hand, on account of self-sow each gene haploid has its rebirth rate: 

N A (t + At) = u A N A (t),N B (t + At) = uj B N B (t) (3) 

Where lo a = 1 + ^-,u B = 1 — If, St is selection gene and At is time gap between border 
generations. From Eq.(l) and Eq.(3) we can get: 

*-">-*> = f^f^}- w 

Considering Eq.(2) and Eq.(4) together we can get the differential equation of x(t) at 
the case of At — ► 0. 

x = P — 7X + fxx(l — x). (5) 

Where f3 = m B ,^f = m A + m B ,[i = In order to simplify the equation, we suppose 
m A + m B = 1. So the simplified gene selection dynamic equation is shown as following 



[14] 

x = (3 — x + [xx(l — x). (6) 

Now if due to some environmental external disturbance the gene selection rate of the 
haploid, it is likely to affect both [3 and fi in the form of additive and multiplicatively 
noises that are connected through a correlation parameter. In other words the external 
fluctuations affect the parameter which fluctuates around a mean value, thus generating 
multiplicative noise and at the same time environmental fluctuations perturbs the dynam- 
ics directly which gives rise to additive noise. So We have 

x = (3 - x + ^x{\ -x) + x(l - x)e(t) + T(t). (7) 
Where e(t),T(t) are Gaussian white noises with the following properties. 

< e (t) >=< T{t) >= 0, (8) 



<e(t)e(t) >=2D6(t-t), (9) 

< T(t)T(t) >= 2aS(t - t), (10) 

< e{t)T(t) >= 2Av A Da<5(t-t'), (11) 

where D,a are the strength of noise e{t) and T(t), respectively. A denotes the degree 
of correlation between noises e(t) and T(t) with < A < 1. 



3 Steady State Analysis of the Model 

We can derive the corresponding Fokker-Planck equation for evolution of Steady Proba- 
bility Distribution Function (SPDF)based on Eq.(7)-Eq.(ll). The equation is shown as 
follow. [15]. 

dP(x,t) dA(x)P(x,t) d 2 B(x)P(x,t) 

dt dx dx 2 { ' 

Where 

A(x) = 2Dx 3 - (fi + 3D)x 2 + (D + fi - 1 - 2\VDa)x + p + X^Da. (13) 

B{x) = D[x{l - x)] 2 + 2Av A Dax(l - x) + a. (14) 



3 



The steady probability distribution of Fokker-Planck equation is given by [14] 



Pst{x) 



B{x) 



exp| 



A(x) 
B{7) 



dx 



(15) 



Where iVo is the normalization constant, using the forms of A(x) and B(x) we have 
the following integral forms of the SPDF of Eq.(12)[16]. 



Pst{x) 



N 



D[x - x 2 )] 2 + 2Av A Da(x - x 2 ) + a 



: exp 



/2, 



\J3 — x + /j,(x — x )]dx 



D{x - x' 2 ) 2 + 2\^fD^(x - x' 2 ) + a 



(16) 

Now we will give a numerical analysis for Eq. (16) and the results will be presented as 
Fig. 1-Fig. 4. 




Figure 1: Plot of P s t(x) (denotes the probability) against x (the ratio of gene A number 
to the total) for different values of the selection gene \i. -D=0.5,a=2.0,A=0.5,/3=0.5 and 
fj,=-l. 0,-0. 5,0.0,0. 5 and 1.0, respectively, (units are relative) 

In the Fig.l, we study the effect of the selection rate /x on SPDF. When the selection 
rate is zero (n = 0), the curve is symmetry. It is said that gene A haploid (x) and gene B 
haploid (1 — x) have the same probability distribution. If the selection rate get a negative 
value the peak of the curve is biased to the left. It is evident that the selection is unfair 
which is propitious to gene B haploid. On the other hand, the selection is propitious to 
gene A in case of fi > 0. In order to study the effect of the noise on the system easily we 
adopt the symmetrical case for the following discussion. 
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Figure 2: Plot of P s t{x) (denotes the probability) against x (the ratio of gene A number 
to the total) for different values of multiplicative noise intensity D. /i=0,a=0,A=0,/3=0.5 
and .0=0.5,2.0 and 8.0, respectively, (units are relative) 

In Fig. (2) we show the effect of multiplicative noise on SPDF. For a small value of 
D the curve shows a single peak (at x = 0.5) region which indicates that environmental 
selection gives the same chance to both gene A haploid and gene B haploid and it is not 
easy for us to select one type haploid from the group. As the value of D increases, the 
single peak region vanishes and it evolves a double peak region. The left peak is near the 
position x = while the right one is near the position x = 1, which shows that environment 
selection gives a absolute big chance to one of the two (gene A or gene B). In this case 
it is easy to select one type haploid from the group since the other type haploid number 
can be neglected with regard to the selected type number. On the other hand, from the 
Fig. 2 we can know that the multiplicative noise can induce the phase transition from one 
peak to double peak. 

The effect of the additive noise on SPDF is shown as Fig.3. For a small value of a the 
curve represents a single peak region which vanishes for larger value of a. The position 
of the peak is weakly affected by the strength of a, however, its height may be affected 
intensively by a. It is said that the additive noise can not separate the single peak while 
it can make the peak disappear as a diffuse term. 

The curve of the Fig.4 shows the effect of correlation parameter on SPDF. When A 
is zero, namely no correlation, the curve gives one peak region whose peak position is at 
x = 0.5 which shows that the selection is fair to both gene A haploid and gene B haploid. 
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Figure 3: P s t(x) (denotes the probability) against x (the ratio of gene A number to 
the total) for different values of additive noise intensity a. /x=0,L>=0,A=0,/3=0.5 and 
a=0.5,1.0,3.0 and 10.0, respectively, (units are relative) 

As the value of A increases the peak at x = 0.5 vanishes gradually and finally it becomes 
a concave peak. At the same time, the height of the curve at the both side becomes 
larger gradually. In others word, the noise correlation make the probability distribution 
centralize on both sides (x = and x = 1). It is evident that the noise correlation make 
it possible to separate one type haploid from a haploid group. 



4 Conclusion and Summary 

In this paper, the steady state properties of the gene selection model is investigated in 
presence of the correlated Gaussian white noise. We study the effect of the white noise on 
SPDF in the symmetrical case. The multiplicative noise can break the single peak state to 
the double peak state (see Fig.2), which indicates that the selection in the case is prepon- 
derant for one of the two ( gene A or gene B). The additive noise play a diffusing role in 
the process, which causes the probabilities distribution to a equal distribution. The noise 
correlation can also break the symmetry of AB, which induces the especial gene haploid 
to be selected. 

From the above, it is found that the noise can change the nature selection from a equal 
probability selection to a differential probability one, which benefits us to select one type 
haploid from the haploid group. The breaking of the symmetry and the especial gene 
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Figure 4: P s t(x) (denotes the probability) against x (the ratio of gene A number to the 
total) for different values of noise correlation parameter A. / u=0.0,a=2.0,L>=0.5,/3=0.5 
and A=0. 0,0. 2, 0.4, 0.7 and 1.0, respectively, (units are relative) 

selection are very important to produce of gene order and biology evolution. 

On the other hand, the stochastic force can induce the phase transition in our system. 
This viewpoint is completely novel as the traditional viewpoint thinks the stochastic force 
disturbs the phase transition. 
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